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FOREWORD 


This report is submitted to the National Aeronautics and 
Space Administration, George C. Marshall Space Flight Center, 
Huntsville, Alabama, in accordance with the requirements of 
NAS8-30608. 



ABSTRACT 

A method of incorporating a detailed solar cell radiation de- 
gradation model into a convenient computational scheme suitable 
for the Solar Electric Propulsion System is outlined. The study 
shows that several existing codes may be applied in sequence to 
solve the problem. 
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1. INTRODUCTION 

Solar cells used to supply power for space missions are sub- 
ject to degradation by energetic radiations encountered in space. The 
dominant radiation sources are electrons and protons in the trapped 
radiation belts, or protons and heavier nuclei due to solar flares. De- 
gradation due to galactic cosmic radiation is usually negligible. The 
effects of nuclear weapon radiation will not be discussed in this report. 
Damage to adhesives and darkening of cover glass materials due to 
electromagnetic radiation are also excluded. The latter effects are 
often slow and linear with exposure. 

Radiation damage to solar cells depends upon many factors. 
This complexity forces most analyses and predictive methods to use 
simplifying assumptions. The inaccuracies thus introduced force the 
designer to use conservative safety margins for mission planning. De- 
gradation of existing solor cell array power supplies in space provide 
checkpoints, but the rapid evaluation of solar cell types and the uncer- 
tainties of extrapolating these data to other trajectories pose severe 
problems. 

Parameters used to assess cell performance include the open 
circuit voltage j Voc; the short circuit current, Isc; and the maximum 
power available, Pmax. These parameters are easier to deal with than 
the voltage - current performance curve. However, these parameters 
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are functions of light spectrum in the undamaged state. Following radia- 
tion exposure, the parameters are functions of particle type, particle 
energy, particle flue nee, light spectrum, time, and temperature -vs 
time profile. Certain other effects such as the angular distribution of 
the radiation and non-uniform damage across the cell due to low energy 
protons complicate the problem. 

The George C. Marshall Space Flight Center recognized the 

need for more accurate degradation predictive methods for the Solar 

Electric Propulion Stage (SEPS). Accordingly, MSFC commissioned 

Horne and Wilkinson to survey the state-of-the-art and recommend an 

approach to improve damage prediction 1 . Horne and Wilkinson conclude 

that a reasonably detailed transport calculation be performed within the 

cell. A microscopic defect damage coefficient may be calculated at 

various points across the cell. The resulting damage profiles should 

then be analyzed with the PN code developed by GGA under Air Force 
2,3 

contract . 

The present effort attempts to integrate the above recom- 
mendations into a computational scheme. Implementation of the scheme 
has not been possible due to limited funds and a delay in delivery of the 
PN code to MSFC. 
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2. COMPUTATIONAL MODEL 

The computational model devised for solar cell radiation 
degradation models requires several steps. Each stage in the calcula- 
tion is outlined below. 

2. 1 RADIATION ENVIRONMENT 

Trapped radiation levels comprizing protons and electrons 

4 

are described by the models of Vette and his co-workers . A trajectory 
code presently operational at MSFG may be used to estimate radiation 
fluences over specified orbits. 

The estimation of solar radiation fluxes is more difficult. It 

5 

is suggested that the work of Burrell be used to determine the number 
of solar events during a mission at various probability levels. It will 
be necessary to assume a spectral shape and normalize it to event mag- 
nitude. For missions within the magnetosphere, a code developed by 
6 

Wright may be used to estimate flux reduction due to magnetic shielding. 

2. 2 PROTON AND ELECTRON TRANSPORT 

It is necessary to find the differential proton and electron 
spectra at the two faces of the solar cell and at a number of mesh planes 
between the faces. A total of three to five calculations should be suffi- 
cient; the spectra at other postions may be obtained by semi -logarithmic 
interpolation for such thin components. 
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For some satellites, the solar cell arrays are attached to 
a reasonably massive body structure. In such cases the radiation com- 
ponent entering the back face may usually be ignored. The SEPS solar 
cell array will probably be unfurled on a thin backing material which 
may be comparable in thickness to the cell itself. Here it will be nec- 
essary to compute radiation levels entering the back side through the 
attachment sheet as well as that entering the front side through the 
cover material. The front and back radiation components will be 
summed. 

. rr 

Horne and Wilkinson suggest the SPARES code system for 
performing the proton and electron transport calculation. These pro- 
grams are suitable for such use and should give accurate results. 
Alternatively, the burden of converting those programs from an IBM 

installation to the Uni vac computers at MSFC may be avoided by using 

8 in 

codes developed by Burrell and Watts or by Hill et al . which are 
operational at MSFC. Should the latter approach be chosen, two code 
modifications will be desirable. First, the treatment of low energy 
protons should be improved by limited parameter fits to range date 
in the low energy regime. Second, the monodirectional proton beam 
technique should be generalized to isotropic incidence. 
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♦ (Note: Recombination Center Curves have been normalized 

to respective damage coefficient curves.) 


Figure 2: Comparison of Energy dependence of Damage Coefficients 
for Different Types of Solar Cells with the Energy De- 
pendence of Recombination Center Introduction Rates 
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to the low electron mass. The solid data points in Figure 2 are de- 
rived from the assumption that all the recombination centers can be 
characterized by a single energy level for a specific material. This 
assumption greatly simplifies the analysis but may lead to some in- 
accuracies in predicting annealing rates. 

2. 4 SOLAR CELL DEGRADATION CALCULATIONS 

The defect density profiles derived from the method out- 
lined in Section 2. 3 identify the radiation-associated variables. Other 
variables such as light transmittance, impurity profiles, diffusion 
length, etc. , must be derived from experimental data and supplied 
to the PN code. This code may then be used to predict gross behavior 
of the solar cell based upon microscopic treatment of the interactions 
Within the cell. 

At the end of the present contract period, the PN code has 
not been delivered to MSFC. This delay has precluded the develop- 
ment of efficient approaches to SEPS mission analysis based on para- 
metric rims. At this point it is not known whether the superposition 
principle will be valid for extrapolating damage calculations to long 
time intervals, particularly for steep damage profiles within the cell. 
The question of annealing during a fluctuating temperature profile is 
also unresolved. A study of these factors must be made in order to 
devise a rapid analysis tool suitable to the SEPS requirements. 
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3. RECOMMENDATIONS 

In order to advance the development of a radiation degrada- 
tion analysis tool suitable for SEPS missions, the following steps are 
recommended, 

1. Increase the accuracy of the proton range function parameters down 
to 0. 2 MeV for the MSFC transport code. The code should be modified 
so that it can sum isotropic proton and electron fluxes incident from 
both sides of the array, fold in the defect cross sections, and output 
the displacement defect profiles. 

2. Conduct a parametric damage study using the PN code in order 
to determine degradation sensitivity to critical parameters. Know- 
ledge of these sensitivities will permit a choice of suitable approxi- 
mations to be made in the analysis method. 


9 



REFERENCES 


1. Horne, W. E. and Wilkinson, M. C. : Research Study on High 
Energy Radiation Effect and Environment - Solar Ceil Degrada- 
tion Methods, D1 80-18475-1, Boeing Aerospace Company Re- 
port, October 24, 1974. 

2. R. E. Leadon and M. L. Vaugh, ’’Short-Pulsed Radiation Effects 
on Dynamic Electronic Components”, Final Report, Contract 
DASA 01-68-C-0123, Prepared for Defense Atomic Support 
Agency by Gulf General Atomic Incorporated, June, 1969. 

3. R. E. Leadon, et al. , ’’Radiation Effects on Dynamic Electronic 
Components, Final Report, Contract DASA 01-69-C-0113, Pre- 
pared for Defense Atomic Support Agency by Gulf Radiation 
Technology, July, 1970. 

4. J. I. Vette, et al. , NASA SP 3024, Vol. 1, 2, . . . , 6 and NS 
SDC 72-06, 72-13, NASA Reports from 1965 to the present. 

5. Burrell, M. O. : The Risk of Solar Proton Events to Space 
Missions, NASA TN D-6379, June 1971. 

6. Burrell, M. O. and Wright, J. J. : The Estimation of Galactic 
Cosmic Ray Penetration and Dose Rates, NASA TN D-6600, March 
1972. 

7. Paul Hahn, ”RSIC Computer Code Collection- -Space Radiation 
Environment Shielding Systems (SPARES), CCC-148, AS-2807, 

Oak Ridge National Laboratory Radiation Shielding Information 
Center, Oak Ridge, Tenn. , 1969. 

8. Burrell, M. O. : The Calculation of Proton Penetration and Dose 
Rates, NASA TM X-53063, 1964. 

9. Watts Jr., J. W. and Burrell, M. O. : Electron and Bre ms - 
strahlung Penetration and Dose Calculation, NASA TN D-6385, 
June 1971. 

10. Hill, C. W. , ; Ritchie, W. B. ; Simpson Jr. , K. M. : Data 
Compilation and Evaluation of Space Shielding Problems, Vol. 

IV - LSVDC4 Program System, Lockheed Report ER-7777, 1967. 


10 



11. J. Baicker, et al. , Applied Physics Letters, 2, 104, 1963. 

12. J. R. Carter, Effect of Electron Energy on Defect Introduction 
in Silicon’*, J. Phys. Chem. Solids, 27, 913, 1966. 


11 



